


use "~/05_vreeland_2008/replication_data.dta", clear
* Change path as needed



*** Table 3: The effect of multiple parties on torture in dictatorships
* Ordinal logit

* Original: Vreeland's table 1, model 1
ologit torture party level g population tradegdp war communist if reg==1, robust
eststo model1, title("Original Model")

* Using 1998 data
ologit torture party level_1998 g_1998 population_1998 tradegdp_1998 war communist if reg==1, robust
eststo model2, title("WDI 1998")

* Using 2004 data
ologit torture party level_2004 g_2004 population_2004 tradegdp_2004 war communist if reg==1, robust
eststo model3, title("WDI 2004")

* Using 2018 data
ologit torture party level_2018 g_2018 population_2018 tradegdp_2018 war communist if reg==1, robust
eststo model4, title("WDI 2018")

* Export
esttab model1 model2 model3 model4 using results.tex, mtitle star(* 0.1 ** 0.05 *** 0.01) b(3) se scalars("ll Log Likelihood") nogaps replace





******** APPENDIX ********

use "~/05_vreeland_2008/replication_data_complete_observations.dta", clear
* Change path as needed

*** Table B.1: The effect of multiple parties on torture in dictatorships
* Ordinal logit
* Reducing all vintages to the common denominator

* Original: Vreeland's table 1, model 1
ologit torture party level g population tradegdp war communist if reg==1, robust
eststo model5, title("Original Model")

* Using 1998 data
ologit torture party level_1998 g_1998 population_1998 tradegdp_1998 war communist if reg==1, robust
eststo model6, title("WDI 1998")

* Using 2004 data
ologit torture party level_2004 g_2004 population_2004 tradegdp_2004 war communist if reg==1, robust
eststo model7, title("WDI 2004")

* Using 2018 data
ologit torture party level_2018 g_2018 population_2018 tradegdp_2018 war communist if reg==1, robust
eststo model8, title("WDI 2018")

* Export
esttab model5 model6 model7 model8 using results.tex, mtitle star(* 0.1 ** 0.05 *** 0.01) b(3) se scalars("ll Log Likelihood") nogaps replace




use "~/05_vreeland_2008/replication_data.dta", clear
* Change path as needed

*** Table B.2: The effect of multiple parties on torture in dictatorships
* Fixed effects logit

* Original: Vreeland's table 1, model 2
clogit common_tort_hath party level g population tradegdp war communist if reg==1, group(aclpcode) 
eststo model9, title("Original Model")

* Using 1998 data
clogit common_tort_hath party level_1998 g_1998 population_1998 tradegdp_1998 war communist if reg==1, group(aclpcode) 
eststo model10, title("WDI 1998")

* Using 2004 data
clogit common_tort_hath party level_2004 g_2004 population_2004 tradegdp_2004 war communist if reg==1, group(aclpcode) 
eststo model11, title("WDI 2004")

* Using 2018 data
clogit common_tort_hath party level_2018 g_2018 population_2018 tradegdp_2018 war communist if reg==1, group(aclpcode) 
eststo model12, title("WDI 2018")

* Export
esttab model9 model10 model11 model12 using results.tex, mtitle star(* 0.1 ** 0.05 *** 0.01) b(3) se scalars("ll Log Likelihood") nogaps replace



*** Table B.3: The effect of multiple parties on torture in dictatorships
* Duration dependence logit

* Generate splines
btscs common_tort_hath  year aclpcode, generate (hath_spline) nspline (3)

* Original: Vreeland's table 1, model 3
logit common_tort_hath party level g population tradegdp war communist hath_spline _spline1 _spline2 _spline3 if reg==1, robust
eststo model13, title("Original Model")

* Using 1998 data
logit common_tort_hath party level_1998 g_1998 population_1998 tradegdp_1998 war communist hath_spline _spline1 _spline2 _spline3 if reg==1, robust
eststo model14, title("WDI 1998")

* Using 2004 data
logit common_tort_hath party level_2004 g_2004 population_2004 tradegdp_2004 war communist hath_spline _spline1 _spline2 _spline3 if reg==1, robust
eststo model15, title("WDI 2004")

* Using 2018 data
logit common_tort_hath party level_2018 g_2018 population_2018 tradegdp_2018 war communist hath_spline _spline1 _spline2 _spline3 if reg==1, robust
eststo model16, title("WDI 2018")

* Export
esttab model13 model14 model15 model16 using results.tex, mtitle star(* 0.1 ** 0.05 *** 0.01) b(3) se scalars("ll Log Likelihood") nogaps replace



*** Table B.4: The effect of multiple political parties on CAT signing 
* Weibull hazard model

* Generate indicator of "failure" (signing the CAT, cats==1)
stset year, id(aclpcode) origin(time max(1984,bornyear)) failure(cats==1)

* Original: Vreeland's table 3, model 2
streg party log_tort communist regscore n_cats muslim level population tradegdp if reg==1, cluster(aclpcode)   dist(weibull)
eststo model17, title("Original Model")

* Using 1998 data
streg party log_tort communist regscore n_cats muslim level_1998 population_1998 tradegdp_1998 if reg==1, cluster(aclpcode)   dist(weibull)
eststo model18, title("WDI 1998")

* Using 2004 data
streg party log_tort communist regscore n_cats muslim level_2004 population_2004 tradegdp_2004 if reg==1, cluster(aclpcode)   dist(weibull)
eststo model19, title("WDI 2004")

* Using 2018 data
streg party log_tort communist regscore n_cats muslim level_2018 population_2018 tradegdp_2018 if reg==1, cluster(aclpcode)   dist(weibull)
eststo model20, title("WDI 2018")


* Generate indicator of "failure" (ratifying the CAT, catr==1)
stset year, id(aclpcode) origin(time max(1984,bornyear)) failure(catr==1)

* Original: Vreeland's table 3, model 4
streg party log_tort communist regscore n_catr muslim level population tradegdp if reg==1, cluster(aclpcode)   dist(weibull)
eststo model21, title("Original Model")

* Using 1998 data
streg party log_tort communist regscore n_catr muslim level_1998 population_1998 tradegdp_1998 if reg==1, cluster(aclpcode)   dist(weibull)
eststo model22, title("WDI 1998")

* Using 2004 data
streg party log_tort communist regscore n_catr muslim level_2004 population_2004 tradegdp_2004 if reg==1, cluster(aclpcode)   dist(weibull)
eststo model23, title("WDI 2004")

* Using 2018 data
streg party log_tort communist regscore n_catr muslim level_2018 population_2018 tradegdp_2018 if reg==1, cluster(aclpcode)   dist(weibull)
eststo model24, title("WDI 2018")

* Export
esttab model17 model18 model19 model20 model21 model22 model23 model24 using results.tex, mtitle star(* 0.1 ** 0.05 *** 0.01) b(3) se eform scalars(obs) nogaps replace

